Fractal Space -Time from Spin- Foams 
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In this paper we perform the calculation of the spectral dimension of spacetime in 4d quantum 
gravity using the Barrett-Crane (BC) spinfoam model. We realize this considering a very simple 
decomposition of the 4d spacetime already used in the graviton propagator calculation and we 
introduce a boundary state which selects a classical geometry on the boundary. We obtain that 
the spectral dimension of the spacetime runs from « 2 to 4 when the energy of a probe scalar field 
decreases from high E < Ep /25 to low energy. The spectral dimension at the Planck scale E ~ Ep 
depends on the areas spectrum used in the calculation. For three different spectra + 1), 

l P (2j + 1) and 1% j we find respectively dimension m 2.31, 2.45 and 2.08. 



Introduction. In past years many approaches to quan- 
tum gravity studied the fractal properties of quantum 
spacetime. In particular in causal dynamical triangula- 
tion (CDT) [l[ and asymptotically safe quantum gravity 
(ASQG) 0, a fractal analysis of spacetime gives a two 
dimensional effective manifold at high energy. In both 
approaches the spectral dimension is D s = 2 at small 
scales and T> s = 4 at large scales. The previous ideas 
have been applied also in the context of noncommutative 
geometry to a quantum sphere and K-Minkowski [3[, in 
causal sets [iKH and in loop quantum gravity Sec 
7} for a summary and new insights. Spectral properties 
have been considered also for the cosmology of a Lifshitz 
universe 0]. Spectral analysis is a useful tool to under- 
stand the effective form of space at small and large scales. 
We believe that the fractal analysis could be also a use- 
ful tool to predict the behaviour of n-point correlation 
functions at small scales Q [l(| and to attack the singu- 
larity problems of general relativity in the full theory of 
quantum gravity [llT|[T^|. 

In this paper we apply to the Barrett-Crane (BC) spin 
foam model [l3j for Riemannian 4d quantum gravity the 
analysis introduced in f(|. The fractal properties of 3d 
spacetime were studied in [Tij ]. The main ingredient is 
the general boundary formalism which provides quan- 
tum amplitudes associated to a finite spacetime region 
for a given boundary 3-geometry [l5][16]. The formal- 
ism is implemented in the same spirit of the calculation 
of LQG graviton propagator [l^l @| E3l E3l C^i [23l E2| - 
For our purposes we will consider a boundary state 
peaked over the boundary geometry of a single 4-simplex. 

The paper is organized as follows. In the first sec- 
tion we define the framework and recall the definition 
of spectral dimension in diffusion processes. In the sec- 
ond section we define the spectral dimension in quantum 
gravity. The analysis is general and not strongly related 
to the specific models. We continue calculating explic- 
itly the spectral dimension for the BC theory using the 
general boundary formalism to define the 4d quantum 
gravity path integral. 

a. Spectral dimension in diffusion processes. The 
following definition of fractal dimension is borrowed from 



the theory of diffusion processes on fractals [23| and eas- 
ily adapted to the quantum gravity context. Let us 
study the Brownian motion of a test particle moving on 
a d-dimensional Riemannian manifold M. with a fixed 
smooth metric g^ v {x). The probability density for the 
particle to diffuse from x' to x during the fictitious time 
(this is just a fictitious time since we are probing the 
spacetime properties, not only the properties of space) 
is the heat-kernel K g (x,x';T), which satisfies the heat 
equation 



d T K g (x,x';T) = A g K g {x,x';T) 
where A„ denotes the covariant Laplacian: 
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A g cf>=—d^g^d v 



(1) 



(2) 



The heat-kernel is a matrix element of the oper- 
ator exp(TA g ), acting on the real Hilbert space 
L 2 (A4, ^/gd d x), between position eigenstates 

K g (x,x';T) = (x'\eMTA g )\x). (3) 

Its trace per unit volume, 



P^T^V- 1 / r\ d x^gJx)K g {x,x-T) 



V^ 1 Tr exp(TA s ) 



(4) 



has the interpretation of an average return probability. 
Here V = J d d x ^fg denotes the total volume. It is well 
known that P g (T) possesses an asymptotic expansion for 
T -> of the form P g (T) = {AirT)- d ' 2 J2n=a A n T n . The 
coefficients A n have a geometric meaning, i.e. A is the 
volume of the manifold and if d = 2 then A\ is propor- 
tional to the Euler characteristic. From the knowledge 
of the function P g one can recover the dimensionality of 
the manifold as the limit for small T of 



m dlnPJT) 
dlnT 



(5) 



If we consider arbitrary fictitious times T, this quantity 
may depend on the scale we are probing. ((5|) is the defi- 
nition of fractal dimension we will use. 
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b. Spectral dimension in quantum gravity. In quan- 
tum gravity it is natural to replace the Laplacian in the 
heat-kernel equation ([T]) with an effective Laplacian as- 
sociated to the expectation value of the operator (A g ). 
It is defined by the path integral 



<A § > = / DgA g e iS ^. 



(6) 



Using the effective Laplacian defined above the heat- 
kernel equation (fTJ) becomes 



d T K(x,x';T) = (A g )K{x,x';T). 



(7) 



Given (A g ) we can define the average return probability 
P(T) in analogy with the classical definition, 



P(T) = y- 1 Trexp(T(A § )) 



(8) 



The spectral dimension of quantum spacetime is defined 
as in $5§. 

c. Including quantum gravity effects. We expect the 
effective Laplacian to run with the probed energy scale k 
when quantum gravity effects become important. Above 
some energy scale the quantum spacetime reveals its hazi- 
ness. So we expect that 



(A«) 



F(k)A k 



(9) 



where A^ Q is the Laplacian at a reference infrared mo- 
mentum fco, and F(k) is a scaling function of the energy 
scale. For instance, in the case the i.r. Laplacian is the 
flat one, we have, expading in Fourier modes: 



P(T)= f d d pe TF ^P 2 



(10) 



We argue that the scales k and p should be physically 
identified, since p represents the energy probed by the 
p— th mode of a scalar field. Then we define an average 
return probability which includes quantum gravity effects 
by 



]d np TF(p)p 2 



P QG (T)= / d d pe 
If, as we will find, F(k) cx k 2 ( z ~ x \ then 

/><,«.< 7.= / dV^cx-^ 

1 2z 



(11) 



(12) 



and the fractal dimension is T> s = d/z. 

d. Spectral dimension in spin foams. In the context 
of spin foam models for quantum gravity we can give a 
precise meaning to the formal path integral ([6]) and define 
the spectral dimension in the general boundary formal- 
ism. We introduce a gaussian state |V> q ) peaked on the 
boundary geometry q = (q,p) defined by the 3-metric 
and the conjugate momentum. We can think the bound- 
ary geometry to be the boundary of a (d— l)-dimensional 
ball. The state is Gaussian and symbolically given by: 



*q(f) 



c 



— (s— q) +ips 



(13) 



The amplitude ([6]) can be defined for a general spin foam 
model in the general boundary framework as 



(W|Ag|* q ) _ J2 Sl , S2 W( Sl )( Sl \A g \ S2 )^(s 2 ) 



E fl WWq( S ) 



(14) 



where W(s) codifies the spin foam dynamics. We will 
consider the BC model with vertex amplitude W(s) pro- 
portional to the lOj-symbol {10j}. 

The fundamental ingredient required to study the frac- 
tal properties of quantum spacetime is the the scaling of 
the metric. To this purpose we consider only one sim- 
plex that could be part of a more complicated simplicial 
decomposition. In other words we suppose spacetime to 
be a 3-ball with a boundary 3-sphere and triangulate the 
sphere in a very fine way (the dual of the triangulation 
is a spin- network). Then we consider another 3-ball and 
its boundary 3-sphere but at a smaller scale. Since we 
are considering the 3-sphere at two different scales all the 
representations labeling the boundary spin-network will 
be rescaled conformally. In our simplification the 3-ball 
is approximated by a single 4-simplex and the boundary 
3-sphere by the ten faces of the simplex. Summarizing: 
in 4d gravity we approximate a 3-ball with a single sim- 
plex and the boundary sphere S 3 with the surface of the 
simplex given by ten triangles. 

We use a conforma! metric defined by = g° 
where <7°j, is a background metric and fl 2 is the conformal 
factor. This is the hypothesis made in Q and [l4j. In the 
quantum theory the conformal factor is replaced by the 
area operator and the metric operator is defined by g^ v — 
Ag° v . To extract the fractal properties of spacetime we 
approximate the metric in the Laplacian with the inverse 
of the area operator in 4d. In fact the Casimir operator 
is related to the area spectrum of a triangular face in 
the simplicial decomposition of space. We consider three 
possible area spectra 



A^ = l P \j{j + l)]i 



(15) 



A' 



where j is the SU(2) representation dual to the triangular 
face. The quantity which is necessary to compute (Ag) 



{g^) = {l/A)g^. 



(16) 



The boundary state in the notation above is 

*j(je) (X e~^ i: -i.«> M ' 1 « (, '' 1 ~ i)w ' 9 ~ J ' )+i * i: - 3 '. (17) 

The 10 spins j e label the boundary triangles of a 4- 
simplex. The Gaussian is peaked over the homoge- 
neous spin configuration, namely over an equilateral 4- 
simplex of size of order y/j. The dihedral angles $ = 
arccos(— 1/4) define the boundary extrinsic geometry. M 
is a 10 x 10 matrix with positive definite real part. We 
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FIG. 1: Plot of the amplitude CEHJ against j for 1 < j < 7. 
The black points represent the amplitude r\ lyV\\j j) and 
the red points represent + 1)]' . 



chose M = 1 ioxio 111 the numerical calculation, but other 
choices do not affect qualitatively the analysis. 
The expectation value (fTl)) reads 



10 

= V J2 II dim(je){10j} , (18) 

where we introduced the following notation for the nor- 
malization, rj^ 1 = (yV\^j). In FigQ]we have plotted the 
amplitude (fT5|) as a function of the SU{2) representa- 
tion j. This is also the quantum gravity scale defined by 
>' l P 3- 

For j > 7 the amplitude (fT5)l is well approximated by 
the classical function 



r,(W\l/Al\^ j 



1 



(19) 



The result (|T9|) can be obtained also analytically in the 
large j limit replacing the amplitude W(j) with the ex- 
ponential of the Regge action. More precisely, since the 
boundary state is peaked on large values of the spins, the 
lOj-symbol in the sum (|18p can be approximated with its 
large spin asymptotic formula 



{10j} ~ A(j) COsSReggcC?) + B(j) 



(20) 



where SRogge is the Regge action for a 4-simplex. One of 
the two exponentials in which the cosine can be decom- 
posed suppresses the sum (|18l) through the presence of a 
rapidly oscillating phase. The term B(j) in (|20|) is dom- 
inant but it is suppressed inside the sum (fT5)) through 
a non trivial mechanism (2~ij |. Expanding SRegge up to 
second order the result (|19p is obtainded by a simple 
gaussian integration. We refer to the literature on the 
graviton propagator for this technique. 

In the range 1 < < 4 we have interpolated the quan- 
tity (fT5|) numerically. Data are fitted with the function 



a l j a i where a rs 0.66, and a ~ 0.73, 
r}(W\Z/A*\9/ ° M) 



J 



0.73 ' 



(21) 



What we learnt from the explicit calculation of (fT5)) can 
be summarized as follows, 



^r— for j > 1 (j > 7), 

tfCi+i)]*£ (22) 

\A & / 1 0.66 - - - 



l 1 



for 1 < j < 7 



The result (|22|) can be used to define an effective Lapla- 
cian (formula (HU) 

a; 

(23) 




l)]l/2 AjO 



for 7 < j < j 
for j « 1 (1 



<i< 



4). 



where we introduced the infrared scale jo 3> 1 ; and, by 
definition, = [j (j + 1)] 1/2 - Actually (J23]) is correct 
also for the area spectra A° and Aj but with different 
values of a, a and the proper functions A^ Q and Aj q . 

We introduce here a physical input to put the mo- 
mentum k in our analysis. If we want to observe the 
spacetime with a microscope of resolution I = Ipj 1 ^ 2 

1/2 

(the infrared length is £q := lpj ) we must use a prob- 
ing particle of momentum p ~ 1/Z. The energy scaling 
property of the Laplacian, i.e. the scaling function F(p), 
can be obtained by replacing: I ~ 1/p, Iq ~ 1/po and 
Zp ~ l/Ep, where po is an infrared energy cutoff and Ep 
is the Planck energy. This gives 



F(p) 



c'p 2a 



+ 1 for p < £| 

for ^r<P<E f 



(24) 



The term +1 is added to facilitate the calculations of the 
spectral dimension. Notice that the scaling function F(p) 
represents also by construction the scaling of the inverse 
metric, (g^} p = F(jp){g^) Vo . 

The parameter z introduced in the second section for 
the area operator A a is 




for P < ff 

for <P<E f 



(25) 



The dependence from T in (jTTJ) determines the fractal di- 
mensionality of spacetime via In the limits T — > oo 
and T — > where we are probing very large and small dis- 
tances, respectively, we obtain the dimensionalities cor- 
responding to the largest and smallest length scales pos- 
sible. These limits are governed by the behaviour of F(p) 
for p — > and p — * oo, respectively. The spectral dimen- 
sion increases from 2 to 4 when the momentum of the 
probe field decreases from p sw Ep/y7 to p w po- A plot 
of the fractal dimension is given in FigJ^l In the regime 
Ep/2 < p < _Ep the spectral dimension is 2? s w 2.31. 
We can repeat the analysis for the area spectra A^ and 
A C y The spectral dimension is the same for j > 5 or 
p < Ep/25 but we can have differences at the Planck 
scale. We summarize the result in the following table, 
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FIG. 2: Picture on the left: plot of the spectral dimension 
f° r i ^ 7 or fc < Ep/y/l as function of the fictitious time 
T. The dimension at close to the Planck scale is 2. We have 
plotted T e [0, 100] and used E p = 100, k = 1. Picture on 
the right: this is an artistic representation of the fractal prop- 
erties of the spacetime. At any scale the spacetime appears 
as a ball (because we consider a compact space time but the 
result is independent from the spacetime topology) and the 
collection of all those balls (one for any scale) originates a 
fractal quantum spacetime. 



Area Spectrum 


z 


Fractal Dimension 




1.73 


2.31 


A* = l P (2j + l) 


1.62 


2.45 


A*, = l 2 P j 


1.92 


2.08 



Correlations. The 2-point correlation function of the 
Brownian motion over a Riemannian d-dimensional man- 
ifold is defined by 



r+oo 

G{x,x') = dTK(x,x';T). (26) 
Jo 

The heat kernel has an approximate expression for x ~ 
x': 

K ^ T ) - jizkm^^i 1 + cu --), (27) 



where a(x,x') is the geodesic distance and with "curv." 
we have indicated the curvature corrections. For x ~ x' 
the result of the integration (|2"6f is 

GM^^M^rg-i). (28) 

We use the spectral dimension we found in the spin foam 
model as an effective dimension at high energies, hence 
in (|28p we make the substitution d — > T> s . At low en- 
ergy p po <C Ep the spectral dimension is T> s = 4 and 
G(x,x') w l/4:Tr 2 a(x, x'). In the energy range Ep/2 < 
k ^ Ep where the spectral dimension is T> s = 2 + e 
the propagator is: G(x, x') w — \og[a{x,x')Tre]/A-K. At 
the Planck scale instead, defining T> s = 2 + (/3 — 
0.31,0.45,0.08 for the three different area spectrum), 
G(x,x') « T{p/2)/(4TT 1+ ^ 2 a(x,x')^ 2 ). 

Conclusions and Discussion. In this paper we com- 
puted the spectral dimension (T> s ) of 4d quantum space- 
time in the Barrett-Crane spinfoam model. We consid- 
ered the simplest decomposition of spacetime and used 
the general boundary formalism to characterize the scal- 
ing properties of the expectation value for the traced 
propagation kernel. Our main observable is a conformal 
metric defined by = il 2 g°, where g°^ v is a back- 
ground metric. In the quantum theory the conformal 
factor is quantized as a function of the area operator. We 
made our analysis with three kind of area spectra, obtain- 
ing different values for the spectral dimension at Planck 
scale (see table) but the same result for E < Ep/2 (see 
Fig®. 

We interpret the results in the following way. At high 
energy the spectral dimension is T> s sa 2 because the man- 
ifold presents holes typical of an atomic structure. 
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